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ABSTRACT 


Infrasonic  Ambient  Ocean  Noise  Spectra  from  Freely  Drifting  Sensors 

by 

Richard  Lee  Culver 
Master  of  Science  in  Oceanography 
University  of  California,  San  Diego,  1985 
Dr.  William  S.  Hodgkiss,  Chairman 

Self-contained,  freely-drifting  Swallow  floats  capable  of  recording  very  low  frequency  (VLF) 
ambient  ocean  noise  are  under  development  at  the  Marine  Physical  Laboratory,  Scripps  Institu¬ 
tion  of  Oceanography,  San  Diego,  California.  The  buoys  are  ballasted  to  neutral  buoyancy  at 
midwater  depth  where  they  record  the  components  of  particle  velocity  from  which  sound  pres¬ 
sure  levels  may  be  derived.  VLF  ambient  ocean  noise  measurements  taken  during  July  1983 
west  of  San  Diego,  California  at  a  depth  of  approximately  1400  meters  compare  favorably  with 
measurements  made  using  other  types  of  sensors.  Recommendations  for  future  work  are  made. 


vii 


Introduction 


Under  Office  of  Naval  Research  sponsorship,  the  Marine  Physical  Laboratory  has 
been  involved  in  the  design,  fabrication  and  testing  of  prototype,  self-contained  Swallow  floats 
which  can  record  very  low  frequency  (VLF)  ambient  ocean  noise  over  extended  periods  of  time. 
The  autonomous  buoys  measure  and  record  the  components  of  particle  velocity  in  the  1-10 
Hz  band.  They  can  be  ballasted  to  neutral  buoyancy  at  a  desired  depth.  The  buoys  generate 
and  receive  high  frequency  acoustic  signals  which  may  be  used  to  determine  their  relative  posi¬ 
tions.  The  deployment  of  several  of  these  units  forms  a  freely  drifting  array  of  sensors.  Indivi¬ 
dual  element  time  series  may  be  combined  coherently  off-line  using  a  beamformer  (e.g.  the 
MPL  dynamic  Beamformer  [  1  ]). 

The  Swallow  float  design  minimizes  self-noise  which  can  limit  accurate  ambient  ocean 
noise  measurements.  The  floats  drift  freely  and  are  not  subject  to  flow  noise  or  cable  strum¬ 
ming.  They  measure  particle  velocity  and  are  therefore  insensitive  to  variations  in  local  pres¬ 
sure. 

MPL  has  conducted  Swallow  float  deployments  annually  since  1982.  The  1982 
deployment  tested  only  the  acoustic  positioning  system.  Between  12  and  14  July  1983,  three 
Swallow  floats  were  deployed  at  a  location  approximately  50  miles  west  of  San  Diego,  Califor¬ 
nia  and  southeast  of  San  Clemente  Island  (32  °N,  118  °W).  Sea  state  during  the  deployment 
varied  between  1  and  2.  Water  depth  was  approximately  1900  meters.  Two  of  the  floats 
deployed  to  a  depth  of  about  1400  meters  and  the  third  deployed  to  the  bottom.  Measure¬ 
ments  recorded  by  float  #5  which  deployed  to  1400  meters  have  been  analyzed  and  results  are 


reported  in  this  thesis. 


I.  The  Infrasonic  Sensor  Buoy  and  Array 


The  Swallow  float  buoy  is  a  17  inch  diameter  glass  float  containing  three  geophones 
as  directional  velocity  transducers,  a  compass  for  detecting  buoy  heading,  an  acoustic  tran¬ 
sponder  for  positioning,  a  solid  state  memory  data  buffer,  a  digital  tape  data  recorder,  and  an 
acoustically  actuated  ballast  release.  General  hardware  configuration  is  shown  is  Figure  1.1. 

While  deployed,  each  buoy  internally  records  compass  heading,  acoustic  positioning 
pulse  arrival  times  and  three  channels  of  geophone  data.  Deployment  time  is  divided  into  90 
second  periods  consisting  of  89  seconds  during  which  data  are  sampled  and  stored  in  the  buffer 
and  1  second  during  which  the  buffer  contents  are  written  to  tape.  No  geophone  data  are 
sampled  during  the  1  second  buffer  dump.  Tape  recorder  capacity  limits  submergence  time  to 
about  60  hours. 

Once  a  buoy  has  surfaced  and  been  retrieved,  its  data  tape  is  read  using  a  computer. 
Data  from  each  90  second  period  are  organized  in  one  record.  Buoy  heading  and  other  header 
information  may  be  extracted  from  each  record  to  determine  how  those  quantities  varied  dur¬ 
ing  deployment.  Positioning  pulse  arrival  times  can  be  used  to  calculate  ranges  to  the  sea  sur¬ 
face  and  bottom,  other  objects  or  other  pulse  emitters.  Geophone  data  provide  information 
about  the  buoy’s  movement  and  can  be  used  to  derive  information  about  the  sound  field  where 
the  buoy  was  deployed. 

When  several  Swallow  floats  are  deployed  to  form  a  randomly  distributed  array,  their 
relative  positions  as  functions  of  time  can  be  estimated  from  positioning  pulse  arrival  times 
recorded  by  the  individual  buoys.  An  absolute  coordinate  system  can  be  established  using 
receivers  or  transmitters  at  known  positions  (e.g.  a  ship  positioned  via  satellite  navigation  or  a 
bottom  tethered  buoy  whose  position  is  accurately  known).  Geophone  data  can  be  coherently 
combined  using  a  beamformer  capable  of  processing  data  from  directional  sensors  which  are 
randomly  and  dynamically  spaced  and  randomly  oriented.  Using  this  sensor  array,  beams  may 
be  formed  in  the  horizontal  and  verticle  directions  to  investigate  the  directionality  and 


statistical  characteristics  of  ambient  ocean  noise. 


Figure  1.1 


II.  The  Acoustic  Positioning  System 


The  Swallow  float  acoustic  positioning  system  operates  at  8  kHz,  has  a  maximum 
range  of  about  15  km  and  has  a  nominal  resolution  of  1.5  meters  based  upon  a  1  kHz  sampling 
rate  at  the  output  of  the  pulse  arrival  detection  circuit.  Data  from  1982  and  1983  Swallow 
float  deployments  demonstrate  the  system’s  capability  to  measure  buoy  depth. 

A.  System  Functioning 

System  operation  is  depicted  is  Figure  2.1.  One  buoy  emits  a  pulse  each  90  second 
period.  The  buoys  take  turns  pinging  until  all  have  transmitted  and  the  sequence  repeats.  All 
buoys  record  pulse  arrival  times  during  every  period.  Although  only  one  pulse  is  transmitted 
during  each  period,  the  buoys  receives  many  pulses  due  to  reflection  and  reverberation. 

The  arrival  time  of  received  pulses  is  estimated  using  a  pulse  detection  circuit  whose 
output  is  sampled  at  a  1  kHz  rate.  A  bit  stream  is  generated  composed  of  one’s  indicating  the 
presence  of  a  pulse  and  zero’s  indicating  that  no  pulse  is  present.  The  bit  stream  is  inspected 
in  blocks  of  8  bits.  If  any  bits  are  found  set  to  one,  the  block  is  recorded  by  the  buoy  as  one 
byte  along  with  its  associated  time  (two  bytes).  Otherwise,  nothing  is  recorded.  Each  buoy 
can  record  up  to  85  3-byte  groups  during  each  1.5  minute  period.  Since  each  buoy  maintains 
its  own  internal  clock,  deciphering  which  buoy  transmitted  during  a  particular  interval  is  not 
difficult. 

B.  System  Measurements 

Acoustic  positioning  system  surface  echo  data  have  been  shown  to  provide  a  means  of 
determining  buoy  depth  as  a  function  of  time  [  2  ].  Figure  2.2  shows  surface  echo  data  from 
the  July  1982  deployment  in  which  one  Swallow  float  was  deployed  to  a  depth  of  about  900 
meters.  The  uppermost,  dark  edge  of  the  band  which  descends  from  the  left  and  levels  off  as  it 
moves  across  the  plot  represents  the  shortest  path  from  the  float  to  the  surface.  Underneath 


this  band,  a  faint  iine  can  be  seen  to  ascend  from  the  left  and  level  off  as  it  moves  across  the 
plot.  This  line  represents  the  arrival  of  pulses  reflected  from  the  sea  bottom. 

Figure  2.3  shows  surface  echo  data  from  the  1983  deployment  in  which  Swallow  float 
#3  was  deployed  to  a  depth  of  about  1400  meters.  Water  depth  was  about  1900  meters.  The 
dark  band  which  ascends  from  the  left  and  levels  off  as  it  moves  across  the  plot  represents  the 
arrival  of  pulses  reflected  by  the  sea  bottom.  The  band  which  descends  from  the  left  and  stops 
at  the  bottom  of  the  plot  represents  the  arrival  of  surface  reflections. 

The  surface  and  bottom  reflections  have  crossed  in  the  1983  data,  as  would  be 
expected  from  buoy’s  close  proximity  to  the  bottom.  The  surface  reflection  curve  does  not  con¬ 
tinue  across  the  plot  because  the  buoy  had  filled  the  portion  of  its  buffer  used  to  record  poise 
arrival  times  with  bottom  and  reverberation  returns  before  the  surface  echos  arrived. 


Figure  2.1 
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HI.  Geophone  Measurements 


As  discussed  earlier,  the  Swallow  float  sensors  measure  water  particle  velocity  using 
three  orthogonally-mounted  velocity  geophones.  Each  geophone  generates  a  signal  propor¬ 
tional  to  the  component  of  the  velocity  vector  parallel  to  its  axis.  The  signal  is  conditioned, 
digitized  and  archived  on  cassette  tape  within  the  submerged  buoy.  Subsequently,  the  tapes 
are  read  with  a  computer  and  the  signals  analyzed.  Geophone  data  from  buoy  #  3,  recorded 
during  the  July  1983  deployment,  were  analyzed  in  this  way  and  the  results  are  presented  in 
this  section. 

A.  Geophone  Signal  Conditioning 

Signal  conditioning  which  occurs  within  the  Swallow  float  is  depicted  in  Figure  3.1. 


Figure  3.1 


The  z  axis  is  positive  in  the  vertical  direction  and  the  y  axis  is  positive  in  the  direc¬ 
tion  of  the  buoy’s  heading.  Geophone  and  anti-aliasing  filter  frequency  responses  are  shown  in 
Figure  3.2.  The  anti-aliasing  filter  has  a  10  Hz  corner  frequency  and  attenuates  signal  by  15 
dB  at  12.5  Hz  and  50  dB  at  16  Hz.  Each  channel  contains  95  dB  of  fixed  gain  plus  a  variable 
(0  to  26  dB)  gain.  Variable  gain  level  is  controlled  by  the  Automatic  Gain  Control  (AGC)  cir¬ 
cuit,  which  optimizes  utilization  of  the  Analog  to  Digital  Converter  (ADC)  dynamic  range.  Its 
functioning  is  explained  later  in  this  section. 

B.  AGC  Level  and  Buoy  Heading 

AGC  level  is  a  general  indicator  of  geophone  signal  level  and  can  therefore  be  looked 
at  as  a  first  step  in  analyzing  geophone  data.  The  same  AGC  level  is  used  for  all  three  chan¬ 
nels  throughout  any  given  period,  i.e.,  AGC  changes  are  made  only  between  periods.  A  period 
consists  of  89  seconds  during  which  data  are  sampled  and  stored  in  the  buffer  and  1  second 
during  which  the  buffer  contents  are  written  to  tape.  No  geophone  data  are  sampled  during 
the  1  second  buffer  dump.  The  AGC  level  is  constant  throughout  a  period  and  is  based  upon 
samples  taken  during  the  previous  period.  The  following  scheme  is  used: 

1.  If  more  than  1%  of  the  samples  taken  during  the  previous  period  were  at  the  max¬ 
imum  positive  or  negative  value,  the  AGC  level  for  the  following  period  is  stepped 
down  by  1/2  dB. 

2.  If  fewer  than  1%  of  the  samples  taken  during  the  previous  period  were  at  the  max¬ 
imum  positive  or  negative  value,  the  AGC  level  for  the  following  period  is  stepped  up 
by  1/2  dB. 

Thus  AGC  level  goes  down  as  the  geophone  signal  level  goes  up  and  vice  versa. 

AGC  level,  compass  heading  and  battery  voltage  recorded  by  Swallow  Float  #  3  dur¬ 
ing  the  1983  deployment  are  shown  in  Figure  3.3.  (Individual  points  are  not  connected  by  lines 
in  order  to  minimize  the  effect  of  a  few  bad  points.)  AGC  level  decreases  to  zero  during  the 
first  1/2  hour  of  deployment  and  remains  there  until  about  8  hours  after  deployment,  when  it 


begins  to  increase.  This  is  explained  as  follows.  Prior  to  deployment,  the  AGC  level  was  ini¬ 
tialized  to  a  value  of  about  13  dB.  Launching  and  travel  to  depth  were  movement  enough  to 
cause  more  than  1%  of  the  ADC  outputs  to  continually  be  at  the  maximum  positive  or  nega¬ 
tive  level.  The  gain  therefore  decreased  until  it  reached  zero,  where  it  remained  until  the 
buoy  reached  the  depth  at  which  it  was  neutrally  buoyant.  There  its  movement  became 
governed  by  movement  of  the  local  water  particle  field.  This  movement  was  of  small  magni¬ 
tude  so  the  gain  stepped  up  at  the  rate  of  1/2  dB  per  period  until  1%  of  the  samples  were 
again  saturated,  at  which  time  the  gain  stepped  down. 

It  can  be  seen  in  Figure  3.3  that  the  gain  oscillates  within  one  rather  narrow  range  of 
values  during  the  period  between  8  and  13  hours  into  the  deployment  and  within  another 
different  range  of  values  between  13  and  18  hours  It  is  not  completely  understood  why  the 
AGC  level  stabilized  at  these  two  different  values  during  deployment.  One  possibility  is  that 
AGC  level  during  the  13  to  18  hour  period  was  controlled  by  movement  of  the  tape  recorder. 
It  is  a  four-track  serpentine  recorder  which  operates  in  the  forward  direction  for  the  first  13 
hours  of  deployment  and  in  the  reverse  direction  for  the  subsequent  13  hours,  and  the  amount 
of  mechanical  energy  which  it  imparts  the  buoy  when  starting  and/or  stopping  may  depend 
upon  its  direction  of  travel. 

The  buoy’s  compass  heading  varied  only  moderately  during  descent  and  was  rela¬ 
tively  stable  between  about  200  and  090  degrees  throughout  its  stay  at  depth.  The  rather 
strong  correlation  between  buoy  heading  and  AGC  level  during  the  13  to  18  hour  period  may 
be  due  to  buoy  rotation  in  combination  with  the  geophone’s  dipole  beam  pattern  modulating  a 
strong  signal. 

C.  Average  Signal  Power 

As  a  second  step  in  analyzing  geophone  data,  blocks  of  128  consecutive  points  from 
each  geophone  time  series  were  squared  and  averaged  to  form  average  power  sequences.  Fig¬ 
ures  3,4  through  3.7  are  plots  of  average  power  for  selected  periods  of  buoy  #  3’s  1983 


deployment.  Each  figure  represents  15  minutes  (10  records)  of  geophone  data.  Correction  for 
AGC  level  has  been  made. 

The  average  power  plots  contain  several  interesting  features.  Spikes  at  multiples  of 
1.5  minutes  occur  at  record  boundaries  in  most  of  the  segments  and  are  thought  to  be  associ¬ 
ated  with  starting  and  stopping  the  tape  recorder.  The  spikes  are  noticeably  larger  in  the  seg¬ 
ment  at  16.25  hours  and  in  the  second  half  of  the  segment  at  13  hours  than  in  the  other  seg¬ 
ments.  From  Figure  3.3  it  can  be  seen  that  the  smaller  spikes  occur  when  the  AGC  is  stabil¬ 
ized  at  the  higher  level,  while  the  larger  spikes  are  present  when  the  AGC  is  at  the  lower  level. 
The  larger  spikes  may  be  responsible  for  reducing  the  AGC  level,  in  that  they  may  contain 
sufficient  energy  to  cause  more  than  1%  of  the  samples  at  the  output  of  the  ADC  to  be  at 
their  maximum  level  and  thereby  step  the  AGC  level  down.  Spike  size  and  AGC  level  are 
correlated  with  the  direction  of  the  tape  recorder’s  operation. 

Other  spikes  are  evident  at  3.3  and  9.3  minutes  in  Figures  3.4,  3.5  and  3.6,  and  at  0.3, 
6.3  and  12.3  minutes  in  Figure  3.7.  These  spikes  are  caused  by  emission  of  the  ranging  ping  by 
buoy  #  3,  which  occurs  at  20  seconds  into  every  fourth  record.  The  power  level  is  zero  in  two 
sections  of  the  segment  at  13  hours,  indicating  that  geophone  data  could  not  be  retrieved  from 
the  corresponding  records.  Plateaus  are  periodically  evident  in  the  segments  at  12  and  13 
hours.  The  raw  time  series  are  seen  to  be  saturated  when  the  average  power  level  is  at  one  of 
these  plateaus.  It  has  been  found  in  some  instances  that  saturated  time  series  were  caused  by 
a  recorder  capstan  slippage  problem  which  left-shifted  data  by  one  or  more  bits.  (A  left-shift 
of  one  bit  is  equivalent  to  multiplication  by  a  factor  of  2.)  It  is  not  known  whether  this  is 
phenomenon  is  the  cause  of  all  such  periods  of  saturation. 

D.  Geophone  Time  Series 

As  a  next  step,  geophone  data  records  corrected  for  AGC  level  were  investigated. 
Time  series  for  41  second  periods  beginning  9  hours  and  13  hours  after  deployment  are  dep¬ 
icted  in  Figures  3.8  and  3.9  respectively. 


X-axis  data  contain  a  strong  0.28  Hz  signal  added  to  the  background  noise.  This  sig¬ 
nal  is  present  to  a  lesser  degree  in  y-axis  data  and  is  not  observed  in  the  vertical  direction.  At 
0.28  Hz,  the  wavelength  is  about  5400  meters  or  nearly  three  times  the  water  depth.  For  such 
low  frequency  sound,  the  ocean  acts  as  a  wave  guide.  [  3  J  Propagation  involves  complicated 
interaction  with  the  ocean  surface  and  bottom  and  is  very  sensitive  to  the  sound  velocity 
profile  and  the  characteristics  of  the  ocean  bottom.  In  general,  a  wave  pattern  is  generated 
consisting  of  standing  waves  in  the  verticle  direction  and  traveling  waves  in  the  transverse 
direction.  Absence  of  the  signal  in  the  verticle  direction  may  indicate  that  the  buoy  was  posi¬ 
tioned  at  a  node  in  the  standing  wave  pattern. 

The  signal’s  frequency  (17  RPM)  is  low  for  shipping  noise.  [  4  ]  It  has  been  observed  in 
many  but  not  all  of  the  time  series  plots  viewed  thus  far  from  this  deployment  and  its  source 
is  not  known. 

Figure  3.10  depicts  three  channels  of  AGC  gain-corrected  geophone  signals  for  a  41 
second  intervals  beginning  13  hours,  33  minutes  into  the  deployment.  A  boundary  between 
two  90  second  periods  (records  #541  and  #542)  occurs  at  about  7  seconds  on  the  plot  axis.  The 
one  second  period  of  data  writing  has  been  omitted.  The  high  level  transient  present  on  all 
channels  during  the  first  3  seconds  of  the  new  record  is  related  to  the  spikes  in  the  average 
power  sequences  and  is  thought  to  be  caused  by  the  tape  recorder.  The  saturated  levels  at  27 
seconds  on  the  plot  axis  (20  seconds  into  the  new  record),  were  caused  by  emission  of  an  acous¬ 
tic  positioning  pulse  by  this  buoy.  They  also  appear  in  the  average  power  time  series. 

E.  Power  Spectral  Estimates 

As  the  final  step  in  the  analysis,  power  spectra  were  estimated  for  the  time  series 
shown  in  Figures  3.8  and  3.9.  AGC  gain-corrected  data  were  Fourier  transformed  and  a  cali¬ 
bration  curve  was  applied.  The  calibration  curve  corrected  for  fixed  channel  gain  and  fre¬ 
quency  dependent  geophone  and  anti-aliasing  filter  responses. 


Power  spectra  were  estimated  using  the  Welch  method  for  incoherently  averaging 
overlapped  periodigrams.  The  1024  point  data  sequences  were  Fast  Fourier  Transformed 
(FFT)  in  128  point  segments  overlapped  by  50%,  and  the  results  incoherently  averaged.  A 
Kaiser-Bessel  window  (a  —  2.5)  was  applied  to  segments  prior  to  doing  the  FFT.  The  following 
algorithm  was  used: 


G(k)  = 


f,-M-K-U 


\  K  — 1 

E  !  *(*)!  2- 

t=0 


(3.1) 


where: 


G(k)  is  the  value  of  the  k-th  bin  of  the  power  spectrum  estimate, 

f$  is  the  sampling  frequency  in  Hz, 

M  is  the  length  of  the  segment, 

K  is  the  number  of  segments, 

„  M-L 

O 

is  the  power  of  the  window, 


u  -  E  w  (n) 
m 

n= 0 


w(n)  is  the  window  function  sequence, 


M—L 

xi(k)  =  S  w(n)  Xf{n)e 

n= 0 


nk 


is  the  value  of  the  k-th  bin  of  the  FFT 


of  the  ith  segment,  and 

x^(n)  is  the  nth  value  of  the  ith  data  segment. 


The  frequency  of  a  bin  center  may  be  found  from  /  = 


k_ 

{ M 


f9  .  Units  of  (3(£)  are  double-side 


power  spectral  density  and  must  be  multiplied  by  two  when  comparing  with  measurements 
typically  reported  in  the  literature.  They  may  converted  to  decibels  (dB)  by  taking 
10  log  G(k ). 

The  channel  calibration  curve  is  shown  in  Figure  3.11.  It  is  composed  of  the  negative 
of  the  geophone  response  curve  shown  in  Figure  3.2,  the  negative  of  the  anti-aliasing  filter 


response  curve  also  shown  in  Figure  3.2,  a  negative  95  dB  to  back  out  the  fixed  gain  in  each 
channel,  and  a  positive  3  dB  to  account  for  negative  frequencies. 

Figures  3.12  and  3.13  depict  power  spectral  estimates  of  the  time  series  shown  in  Fig¬ 
ures  3.8  and  3.9.  Their  significance  and  relationship  to  sound  pressure  level  is  examined  in  the 


next  section. 
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IV.  Discussion  of  Geophone  Measurements 


Measurements  of  time-varying  particle  velocity  in  the  ocean  are  of  interest  because 
they  can  be  used  to  derive  information  about  the  local  sound  field.  Sound  propagates  through 
a  fluid  by  the  movement  of  fluid  particles,  and  because  the  fluid  is  compressible,  the  particle 
movement  causes  changes  in  pressure.  Using  the  relationship  between  particle  movement  and 
pressure  changes,  certain  characteristics  of  the  sound  wave  can  be  determined  from  a 
knowledge  of  the  particles’  movements. 


A.  Time-Varying  Pressure  and  Particle  Velocity 

The  relationship  between  particle  velocity  and  pressure  changes  can  be  determined 
from  the  physics  which  describe  sound  propagation  in  a  fluid.  The  fundamental  equation 
which  describes  wave  propagation,  the  wave  equation,  may  be  derived  from: 

•  the  equation  of  continuity  (law  of  the  conservation  of  mass), 

•  the  equations  of  motion  (Newton’s  second  law),  and 

•  the  equation  of  state  (relates  pressure  changes  to  density  changes).  [  5  ] 

The  equation  of  continuity  may  be  expressed  as 


da 

dt 


du~  duv  duz 

— -  +  — “  + — L 
dx  dy  dz 


(4.1) 


where: 


a  = ' 


P—Po 

Po 


,  the  fractional  change  in  density, 


/?0  =  the  constant  equilibrium  density, 
p  =  the  actual  density,  and 

uz)uya.nduz  are  the  components  of  fluid  velocity  along  the  x,  y  and  z  axes,  respec¬ 
tively. 
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The  equations  of  motion  may  be  written  as 


fx  =P- 


du * 


dt 


,  ,a”> 

fy  p  dt 


and  fz  =  p- 


duT 


dt 


(4.2) 


where  fx,  fy  and  f z  are  the  force  components  per  unit  volume. 

The  equation  of  state  tells  us  that  the  fractional  change  in  pressure,  denoted  by  p,  is 
directly  proportional  to  a,  the  fractional  change  in  density.  The  constant  of  proportionality  is 
/c,  the  bulk  modulus,  so  that  p  —  kg.  We  also  idealize  the  fluid  as  perfectly  non-viscous, 
requiring  that  it  be  unable  to  support  shear  stresses.  This  assumption  provides  us  with 


fx  = 


dp 
dx  ’ 


and  fz  =  - 


dp 

dz 


(4.3) 


Assuming  that  a  <1,  the  above  equations  can  be  combined  to  form  the  wave  equa¬ 
tion,  which  can  be  written  either  in  terms  of  excess  pressure  or  fluid  velocity 


d2p _ k_ 

dt 2  “  P0 


<Puz  K  d2«s  d2«s  d2tts 

dt 2  Po  dx 2  dy2  dz2  j 


(4.4a) 


(4.4b) 


Substituting  c2  =  —  ,  where  c  is  the  speed  of  sound  in  the  fluid,  yields 
Po 


d^2  _ 

dt2  ~ 


&JL  . 

dx2  dy 2  dz2 1 


and 


d2ux  2  d2ux  d2ux  d2ux 
dt2  ^  dx2  dy 2  dz 2 


(4.5a) 


(4.5b) 


and  corresponding  equations  can  be  written  for  uy  and  uz. 

The  wave  equations  describe  how  pressure  and  velocity  change  as  functions  of  time 
and  position.  Solutions  for  a  plane  wave  traveling  in  an  arbitrary  direction  may  be  written  in 


the  form 
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P+  —  A  F+(ojt  —  KT*r*)  —  A  F+(ut  —  kxx  —  kyy  —  kzz ),  (4.6a) 

P_  — A  F„(gj£  H-F-r")  =A  F+(ut  +  kxz  +  kyy  +kzz)  and  (4.6b) 

P  =  P+  +  P_  =  C  F  (ut  ±  kxx  ±  kyy  ±  kzz ),  (4.6c) 

where  A,  B  and  C  are  constants,  F,  F+  and  F_  are  arbitrary  functions  determined  by  the 

boundary  conditions;  ITis  a  vector  whose  magnitude  is  the  wave  number  (— )  and  which  points 

c 

in  the  direction  of  propagation;  kX)  ky  and  k2  are  the  components  of  IT  along  the  x,  y  and  z 
axes,  respectively;  and  F*is  a  position  vector  which  gives  the  location  of  the  point  ( x,y,z )  with 
respect  to  the  origin,  f  6  ^  Setting 

q  =ut  ±  kxx  ±  kvy  ±  k„z,  (4.7) 

we  can  solve  for 
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(4.8a) 


(4.8b) 


(4.8c) 


(4.8d) 


Substituting  these  relations  into  equation  (4.5)  and  dividing  out  the  common  factor  shows  that 
P,  P+  and  P_  satisfy  the  wave  equation  provided 


kx  +k%  +  kz  — 


r\ 


(4.9) 


To  see  how  pressure  changes  as  a  function  of  the  component  velocity,  the  equations  of  con¬ 
tinuity  and  state  may  be  combined  to  get 
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dt 
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(4.10) 
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Since  ux ,  uy  and  uz  each  satisfy  the  wave  equation  separately,  their  solutions  are  of  the  form 

ux  —  Ax  Fx  (c ot  ±  r*),  (4.11a) 

uy  =Ay  Fy  ( wt  ±  IT'  r*)  and  (4.11b) 

uz  =AZ  Fz  (ojt  ±  F*  r),  (4.11c) 

where  the  A’s  are  constants  and  the  F’s  are  arbitrary  functions  determined  by  the  boundary 
conditions.  With  q  —ut  ±  F*  r*9  we  can  write 
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so  that 
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Substituting  these  relations  into  equation  (4.10)  yields 


(4.13a) 


(4.13b, c) 


dp  _ 
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dt 


dt 


dt 


(4.14) 


which  can  now  be  integrated  to  produce  a  solution  for  the  pressure  in  terms  of  the  component 
velocities, 


(415) 


Substituting  k  =  pQ  c2  and  j  Fj  =  —  produces 
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where 


1  JT1 

i  K  i 


!  rs 


and 


1  JT1 

i  *  i 


(4.17) 


are  called  the  direction  cosines  of  fc*  with  respect  to  the  x,y  and  z  axes. 

Figure  4.1  illustrates  the  special  case  where  the  wave  is  propagating  along  the  x  axis. 
Since  uy  =  uz  =0  and  kx  =  [  F“J  ,  pressure  in  terms  of  particle  velocity  is  p  =  ±  pQcux  . 


Figure  4.1 


Figure  4.2  illustrates  a  second  special  case  in  which  the  wave  is  propagating  along  a  direction 
perpendicular  to  the  z  axis.  For  this  case  uz  =0,  kx  =  |  cos  0  and  ky  —  |  sin  9,  so  that 
the  pressure  is 


± 


=  Po  c( 
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Figure  4.2 


We  model  the  ambient  ocean  noise  field  as  resulting  from  the  arrival  of  an  infinite 
number  of  plane  waves.  At  a  particular  point  in  space  and  time,  the  pressure  resulting  from 
the  ith  plane  wave  is 


Pi  ”  Po  c 
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(4.19) 


and  the  total  pressure  at  a  point  is  the  sum  of  the  pressures  resulting  from  each  wave 
separately, 


p  =  £  Pi  =  Poc  £ 

*=4  *-4 
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(4.20) 


which  says  that  the  velocity  components  for  each  wave  must  be  weighted  with  their  direction 
cosines  and  summed. 

Consistent  with  our  model  of  the  ambient  ocean  noise  field,  the  output  of  the  Swallow 
floats,  ux(n),  uy(n)  and  uz(n),  where  n  is  the  discrete  time  index,  are 

oo 

«*(«)  =  £  «*.(«)>  (4.21a) 
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«*(«)  =  E  «»(«)•  (4.21c) 

t"  5=4 

The  components  of  particle  velocity  in  the  x  direction  due  to  all  waves  have  been  summed. 
Likewise  for  the  particle  velocity  components  in  the  y  and  2  directions.  The  components  asso¬ 
ciated  with  each  wave  cannot  be  extracted  from  the  sum.  Also,  the  angle  of  arrival  of  each 
wave  is  unknown.  We  are  therefore  unable  to  calculate  the  pressure  time  series  p(n)  from  the 
Swallow  float  outputs  tt,(n),  uy{n)  and  uz{n).  Fortunately,  under  certain  assumptions,  we  can 
calculate  the  pressure  spectrum  levels  from  the  component  velocity  spectrum  levels.  This  is 
discussed  in  the  following  section. 


B.  Mean  Square  Power  and  Particle  Velocity 

Using  our  model  of  ambient  ocean  noise  as  resulting  from  the  arrival  of  an  infinite 
number  of  independent  plane  waves,  the  mean  square  pressure  is 
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(4.22) 


Because  the  waves  are  independent,  the  mean  square  total  pressure  is  equal  to  the  sum  of  the 
mean  square  pressures  resulting  from  the  individual  plane  waves.  The  mean  square  pressure 
associated  with  one  plane  wave,  in  terms  of  the  component  velocities  and  using  equation 
(4.20),  is 
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Expressing  the  component  velocities  in  terms  of  the  velocity  vector  of  the  ith  plane  wave  «7* 
and  the  direction  cosines  yields 
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(4.24a, b) 
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Substituting  these  results  into  (4.23)  and  using 


k2.  -U  k 2  A-  k2  —  1  fc? 1  2 
Kxt  *  Kyi  *  Kzt  I  Kt  I 

yields 

e[p?}  =  (/v)2  E  !  “T!  2  =  {p0c)2e  [«£•+«£  + «i<  ]• 
Thus  we  can  write  the  mean  square  pressure  as 


r 

1 

OO  /  \ 

E[i 

2j  =M2E 

E  «£•+«£  +  “«•] 

»-4  V  \ 

(4.25) 


(4.26) 


(4.27) 


Again  invoking  the  independence  of  the  individual  plane  waves,  we  can  equate  the  expected 
value  of  the  sum  of  the  squared  velocity  components  with  the  expected  value  of  the  square  of 
the  sum  of  the  velocity  components,  or 

(4.28) 

This  is  also  true  for  uy  and  uz.  Substituting  these  relations  into  the  previous  equation  yields 

e[p2]  =  (p0c?E[v?  +  «y2  +  «/]•  (4.29) 


oo 

CO 

2  r  i 

E 

E  u% 

E  ux< 
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i=l 

1=4 

The  equation  above  directly  relates  the  mean  square  pressure  to  the  mean  square  of  the 
summed  velocity  components.  It  is  based  upon  our  model  of  the  ambient  ocean  noise  field  as 
resulting  from  the  arrival  of  an  infinite  number  of  independent  plane  waves. 

The  expected  value  of  the  pressure  squared  is  the  average  pressure  power.  On  a  per 
Hertz  basis,  this  is  the  power  spectral  estimate  [7].  Using  the  notation  defined  in  Section  III, 
equation  (4.29)  can  be  written  in  the  form 

Gp(k)  =  (p0cf  [(?.,(*)  +  G,y(k)  +  G,,(*)],  (4.30) 


where  k  is  the  bin  index  of  the  power  spectral  estimate.  Taking  10  log  of  both  sides  yields 
pressure  spectra  in  dB  as  a  function  of  the  velocity  spectra 
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10  logGp(&)  =20  log(p0c)  +10  log  f  Gu  ( k ) 


(4.31) 


For  seawater, 


f  \ 

t  c  v  ^  -i  ^9 

r  \ 

Pa 

10  6nPa 

l.o  X  1U 

mz  •  sec  j 

kg/m  •  sec2  j 

Pa 

(4.32) 


=  1.5XlO*-4^ 

m /sec 


so  that 


10  log Gp(k)  =  243.5  dB  + 10  log  [<?„,(*)  +  G,y{k)  +  G„f(*)]. 


(4.33) 


C.  Swallow  Float  Pressure  Spectra 

Figures  4.3  and  4.4  contain  the  velocity  spectra  shown  earlier  in  Figures  3.12  and 
3.13,  but  with  243.5  dB  added.  In  addition,  the  bottom  plot  in  each  figure  is  the  power  sum  of 
the  three  velocity  spectra  above  it,  which  is  10  log  Gp(£)  in  the  previous  equation.  The  bot¬ 
tom  graphs  in  each  figure  therefore  represent  sound  pressure  levels  derived  from  the  Swallow 
float  particle  velocity  measurements  and  based  upon  the  aforementioned  assumptions  about 
the  composition  of  the  ambient  ocean  noise  field. 

G.  M.  Wenz  [  8  ]  published  six  ambient  ocean  noise  spectra  which  encompass  the  1  to 
10  Hertz  range  and  which  were  acquired  using  omni-directional  hydrophones.  He  did  not  iden¬ 
tify  the  areas  or  depths  where  the  measurements  were  made.  R  H.  Nichols  [  9  ]  has  published 
data  acquired  with  an  omni-directional  hydrophone  bottom-mounted  at  1200  meters  near 
Eleuthera  Island  between  June  and  August  1980.  His  hydrophones  were  enclosed  in  large, 
streamlined  housings  shaped  like  pocket-watch  cases  and  designed  to  shield  the  transducer 
from  self-noise  due  to  current  flow.  Nichols  shows  that  his  data  are  at  or  slightly  below  the 
lower  limit  of  the  envelope  of  the  six  Wenz  spectra.  The  Nichols  data  are  plotted  in  the  bot¬ 
tom  graphs  of  Figures  4.3  and  4.4  for  general  comparison  with  Swallow  float  data. 


Swallow  float-derived  spectrum  levels  are  very  close  to  Nichols’s  measurements. 
Between  1  and  8  Hz,  Swallow  float  levels  display  a  10  -  12  dB  per  octave  downward  slope. 
Levels  rise  sharply  below  1  Hz  and  flatten  above  8  Hz.  These  are  the  same  general  charac¬ 
teristics  which  Wenz  and  Nichols  ascribe  to  ambient  ocean  noise  levels  in  this  band. 

Swallow  float  spectra  appear  rather  smooth  below  below  6  Hz  and  are  somewhat 
bumpy  above.  Ambient  ocean  noise  sources  in  the  1  to  10  Hz  region  are  thought  to  be  pri¬ 
marily  turbulent-pressure  fluctuations.  Above  about  5  Hz,  distant  shipping  noise  may  be  a 
significant  source  [4,8,9].  The  smooth  portion  of  the  Swallow  float  spectra  below  6  Hz  is  there¬ 
fore  most  likely  dominated  by  turbulent-pressure  fluctuations,  while  the  bumpy  features  above 
6  Hz  are  probably  due  to  shipping  noise. 

Turbulence  at  the  ocean  bottom  due  to  current  flow  past  the  ocean-bottom  boundary 
probably  exceeds  that  found  in  the  mid-water  column.  Thus  we  might  expect  spectrum  levels 
measured  by  a  mid-water  column  sensor  such  as  the  Swallow  float  to  be  somewhat  lower  than 
those  measured  at  the  bottom.  Between  5  and  10  Hz,  where  distant  shipping  noise  sources 
begin  to  contribute,  the  difference  in  sensor  location  might  be  expected  to  have  less  effect. 
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V.  Future  Work 

Work  in  the  near  future  will  focus  on  developing  a  capability  to  localize  array  ele¬ 
ments  using  interelement  time  delays.  Data  are  noisy  and  in  some  cases  missing.  One  aspect 
of  the  problem  which  requires  careful  attention  is  how  to  smooth  raw  data  prior  to  using  it  in 
any  localization  calculations.  An  exploratory  look  at  this  problem  is  discussed  in  [  10  ]. 

Subsequently,  element  acoustic  data  may  be  combined  using  a  beamformer  and  the 
directionality  characteristics  of  ambient  ocean  noise  assessed. 
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